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The axial charge qa is calculated in the e regime of Heavy Baryon Chiral Perturbation Theory 
to order e^ . To perform this calculation, we develop a technique to compute baryon properties in 
the e regime of Chiral Perturbation Theory. This technique includes contributions from pion zero 
momentum modes and can be used at arbitrary order, diagram by diagram, in the e regime to 
calculate any matrix element. Also, a calculation of qa in the e' regime to order e'^ is performed. 
A discussion of the domain of applicability for both the t and t regimes is also included. 



I. INTRODUCTION 

As computational speed has steadily advanced in recent years, lattice QCD has been increasingly utilized as the 
only known method of computing observables directly from QCD. However, even the most advanced calculations to 
date use moderate lattice volumes and unphysically large pion masses. State of the art calculations are performed 
with the spatial length of the lattice being on the order of a few Fermis, while the temporal extent is typically larger 
by at most a factor of three. These calculations are also done with pion masses that are at least a factor of two greater 
than the physical pion mass. In studying low energy observables using lattice QCD it is therefore important to fully 
understand how the finite volume of the lattice and the large pion mass affect the result. To study the low energy 
dynamics of QCD and finite volume effects thereof Chiral Perturbation Theory (ChPT) in finite volume can be used. 

When studying baryon properties the low energy theory used is Heavy Baryon Chiral Perturbation Theory 
(HBChPT) y. For infinite volume HBChPT, the small expansion parameters are Aqcd/^b, p/^x' ^^"^ rn^^jK^ 
where -p is the typical momentum and A^ is the chiral symmetry breaking scale typically of order 1 GeV[2]. This 
counting scheme also holds in finite volume when mT^hjl-K » 1[3]. However, a problem exists for small quark masses 
(m7ri/27r << 1), as the zero momentum mode pion propagator goes as l/rn^V, where V is the spacetime volume 4|. 
Thus as computational power continues to increase and lattice calculations are performed at ever lower quark masses, 
the contribution from zero momentum mode pions will become ever larger. To account for this in the HBChPT 
calculation the counting scheme must be changed to enhance the order at which zero mode diagrams contribute. 

To remedy this problem, the construction of a new power counting scheme is needed such that the zero momentum 
modes are enhanced and integrated over exactly while the non-zero momentum modes can be treated perturbatively[4, 
y, y, y, y, y, [31. This leads to the formulation of the e regime. In the e regime, if L and (3 are, respectively, the 
spatial and temporal extent of the box one is working in, e ^ 27r/A^L ^ 27r/A^/3 and TOtt/A^. This counting 
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takes into account that the zero momentum pion contributions have become nonperturbative. The integration over 

. . n 

these zero modes is achieved by utiHzing the method of coUective variables ,4| and methods for including baryons in 
this framework have been proposed previously 

Another important regime for working in small volumes is the e' regime 1^ • This regime is characterized by a 



highly asymmetric hyperbox with spatial dimension L and a large temporal dimension. Due to the zero mode pion 
propagator mentioned above (which goes as as one approaches the e regime from the p regime the spatial 

and temporal zero modes will become enhanced relative to the nonzero modes. While they do not yet need to be 
treated nonperturbatively, their counting does need to be enhanced relative to the nonzero modes. To account for 
this, the e' regime possesses similar counting rules to the e regime: e' ~ 27r/A^L and e'^ ~ m7r/A^, but separate 
rules for the temporal counting. For nonzero momentum modes the time direction counts according to e' ^ 27r/A^/?. 
However, for pion loops containing spatial zero modes this counting changes to e'^ ^ 27r/A^/3. In this regime, zero 
modes remain pcrturbative. For the specific combinations of the lattice size and quark mass used in this paper, the e 
and e' regimes are the correct counting schemes to use. 

The axial- vector current has been previously studied in the p regime in infinite volume . 
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as well as in finite volume 



16l | and provides an important cornerstone of lattice QCD efforts to understand baryon 
24 1 . This work focuses on understanding the nonperturbative aspects of the e regime and 
computing the axial-vector current matrix element between two nucleon states of equal momenta for lattice sizes that 
require the use of the e and e' regimes of HBChPT. By better understanding how the volume of the lattice and choice 
of quark mass determines which power counting scheme is needed, one can more accurately extrapolate the infinite 
volume, physical pion mass results for gA from the lattice data. 



II. HEAVY BARYON CHIRAL PERTURBATION THEORY 



At low energy one can use HBChPT to describe the dynamics of nucleons and pions. The Lagrangian that is 
consistent with spontaneously broken SU{2)l ® SU{2)ji is at leading order [ll [25!: 



f2 f2 
Co = Niv ■ VN - T^iv ■ VT^' + AT^T'' + — Tr[a^E^(9'^I]] + A— Tr[m,I]^ + h.c] 



-^abc,L^ Ad 



(1) 



with the velocity dependent nucleon fields N (for brevity in this paper the usual subscript v and the integral over the 
velocity have been dropped from the nucleon fields), the Rarita-Schwinger fields T'^ describing the A-resonances, and 



the definitions: 

2 f2iM 
E = ^ = exp ( — ^ 

= a^ + v"'^. (2) 

The pion fields are encapsulated in the matrix M . The Rarita-Schwinger fields are rank-3 tensors such that: 

ylll ^ ^++^ 2^112 ^ 2^122 ^ _1^^0^ y222 ^ ^- 

while the nucleons are simply an SU(2) vector given by 

where the constant / = 132 MeV and the matrix rrig is the quark mass matrix. The three couplings given in eqn. 
© (ffA^ 9AN, and qaa) are the infinite volume, chiral limit couplings between baryons and pions. Also, is the 
covariant spin vector and is the heavy baryon four velocity with v'^ ~ 1 (typically = (1,0)). 

In our calculation higher order Lagrangian terms will also become important. The relevant next to leading order 
Lagrangian terms are: 

/-V^-(vV)^ \ /-„2?2_(„.2?)2 

( o'^'^'\ 

+\CiNTi[mqY.'^ + h.c]N + 47V C2 - AIN + ACj.NA^'A^.N 

\ sniB I 

— > N N ~ r„ + ACi^Tr[m„I]^ + h.c.]N + AN \ C2 - AIN + AC^NA^AuN (5) 

2mB 2mB I Sms / 

where the chirally covariant derivatives V are converted to normal derivatives d as the contributions from will not 
be important at the order considered. The coefficients in front of these terms can be determined by reparametrization 
invariance 26] or by matching to the relativistic theory in the path integral 



27j . The last three terms are not important 



for calculating the axial matrix element to order or e'^ , but are important for the nuclcon mass. The first two terms 
in eqn. ([5|) are simple higher order extensions of the kinetic energy operators given in the lowest order Lagrangian. 
Including these terms modifies our baryon and decuplet propagators. It should also be noted that these two order 
l/niB terms are not the only l/rriB operators, but they are the only ones that will enter in the evaluation of the axial 
matrix element to the order we work in e and e'. Defining t^^ = |(C^t°^ + £,t°-£,^) and = ^(^V°^ — S,t°'£,^), the 
axial- vector current is: 



where the elhpses indicate higher order terms. The meson part of the current has been omitted as it is not necessary 
for what follows. The form of the higher order current is easily derivable by beginning with the fully relativistic 
expression for the current, and then substituting in the appropriate expansion for the relativistic field in terms of the 
heavy field to the desired order in 1/niB- Upon doing this, and working to 0{:^), the higher order current is 
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where ]p± = p— ^ {v ■ D), — S'^t^^, and the superscript (1) denotes l/mg suppressed contribution to the current. 
Working to lowest order in the above current amounts to the replacement of p by ^. At tree level, this contribution 
vanishes as it is proportional to the spatial on-shcU momenta which is taken to zero at the end. Similarly, including 
this operator within a loop diagram yields zero. This can be understood simply since both the pion propagator and the 
two pion-baryon vertices are even with respect to the spatial loop momenta whereas the higher order current supplies 
an odd power of the spatial loop momenta. These higher order terms will not be necessary for this calculation. 
In finite volume, the fully relativistic nucleon field must satisfy the conditions llj : 



ijiP,f) = -^(0,r) (8) 

where /3 is the temporal extent of the volume. Choosing the rest frame of the nucleon — (1,0,0,0), this then 
implies that the heavy nucleon field satisfies the condition: 

Ar(/3,f) = -e'^"^iV(0,r). (9) 
This means that the heavy nucleon field has the Fourier decomposition 



N{T,r) = ^exp 



no 



'n{2no + 1) . 
'I 1 h iniB I r 



N{0,r}. (10) 



A similar condition exists for the decuplet fields. The presence of the baryon mass in the above equation will become 
important when looking at the finite time corrections to Feynman graphs. 

III. e AND e' EXPANSION REGIMES OF VALIDITY 

For a nucleon in a finite volume there are four different regimes that are dependent on the choices of the quark mass 
(and by extension the pion mass), the spatial extent of the volume, and the temporal extent of the volume. These 
regimes define particular power countings appropriate for the physics and are variously called e, 6, and e' 



'3,4, 12, 3^ 



regimes, along with the standard p regime. A simple illustration of the boundaries between them can be provided 
by examining how the different modes contained in the simple pion loop in fig. ([1]) count. This definition of the 
boundaries is not unique but serves as a useful guide. Concentrating on the e, S, and e' regimes, each of these three 
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FIG. 1; One Loop Contribution 
expansion parameters are similarly defined by the relation: 

2tt 



A useful relation to separate these regimes from the p regime is the condition 



27r 



< 1 



(11) 



(12) 



which is a quantity that counts as e or e', and allows for a perturbative expansion in m^L. 
Examining the different modes of fig. ([T]) yields: 
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(90,9). 



(13) 



where in finite volume one has the definitions \q\^ — (27r)^(n^ + + n\)/L'^ and — rv^ / 0^ . The spatial and 
temporal zero modes are then defined by n = and tiq — Q respectively. To separate the three regimes (e, 5, and e') 
one must look at the boundaries of where both the spatial and the temporal zero mode enter at order one. 

Also of interest is the quantity AL/27r. It should be noted that strictly speaking AL >> 1 because A ~ Aqcd ~ fir 
and in order for hadronic physics to be contained within the volume we work, f-^L » 1. Calculationally however, 
is numerically the same as in the region of the e and e' regimes that are of primary interest to lattice QCD 
calculations, and this allows for a perturbative expansion of AL/27r. 

In the e regime it is required that the zero mode — (0, 0) counts as order one, while the 5 regime also requires that 
the spatial zero modes = (go, 0) count as order one. In the e' regime all zero mode contributions are perturbative. 
The boundaries of the three regimes in L-rriT^ space for different values of the ratio (3/ L are plotted in fig. ([2]). However, 
it must be emphasized that these boundaries are necessarily very poorly defined in their placement, and should be 
interpreted more as midpoints in the smooth transition between regimes than as hard lines demarcating each region. 

Within each of the regimes, how the pion mass, the decuplet mass splitting, and zero mode (3 count will vary 
depending on how the ratios rriT^/K^, /S./A^^ and 27r/A^/3 compare to the ratio 2'it/K^L. This will in turn affect the 



order at which certain graphs contribute. In this work to^/A^ g2 ^ g/2 A/A^ 



e' . However, it should 



be noted that as one moves toward lattices with smaller spatial dimension and smaller pion mass both m-^/ K^- and 
A/A^ will numerically become smaller relative to 2t: / K~^L and so will count at higher order. This will require a 
reassessment of the order at which the graphs contribute. 
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FIG. 2: Different counting regimes in L-m^ space. The horizontal dashed hne gives the value of 2n/A above which the 
expansion taken in AL is no longer valid. The dashed lines demarcating the boundaries between the different regions are there 
to remind the reader that these are poorly defined boundaries. 

Additionally, the order at which contributes will have nontrivial consequences for how 1/P is counted within 
the e' regime. In this work the zero mode counting of 1/(3 in the e' regime is determined by m^, as this is where the 
zero mode of the tadpole diagram of fig. ([1]) becomes large [l^. However, in assigning 1/(3 ^ the tadpole zero 
mode counts not as \/(3L^m'^ but rather as I/L^ijIt^. This difference determines a boundary within the e' regime 
above which the counting of ^ for a zero mode diagram is legitimate, and below which it is not. This line 
does not demarcate the boundary between e and e', but it does determine the regions where ~ counting is 
appropriate. 

IV. EVALUATION OF THE ZERO MODES FOR THE e REGIME 



For the e regime, the zero modes, = (0,0), must be treated nonperturbatively due to the l/m^V dependance of 
the pion zero mode propagator. Starting with a purely mesonic theory the Euclidean partition function has the form 



Z = j [VT] exp - j d*x CiT,{x)) 



(14) 



The group element T,{x) parameterizes the meson fields. The zero and non-zero modes are described by the Fourier 
components of the pion fields, (7,", in finite volume: 

= exp (^^^^^^M:^^ , Ar(a;)-X^<Z>„(a;), (15) 

where n is a four- vector with integer entries, r° is one-half the Pauli matrices ct", and the u„(a;)'s are plane waves. 
Following the procedure of ref. \^ , a change of variables is made such that the zero modes are separated from the 
non-zero modes: 

E ^ U±U, (16) 



U = exp(^-i^j, (17) 
- , , / 2iM(x) • t\ 

E(a;) = expl 7^)' (18) 

M^{x) = J2 Pn,.^",ix). (19) 

The U^s are spacetime independent SU(2) matrices containing the zero mode contributions while Yl{x) contains the 
non-zero modes with new Fourier coefficients Note the pj^'s are related non-linearly to the g°'s. Using the Baker- 
Campbell-Hausdorff formula, all the exponentials in eqn (|16p can be combined. The resulting exponential has the 
form: 

UtU = exp(^(l)''T'' + ^M''{x)t'' + ^[2irT'',2iM\x)T''] + . 
\ J J 




'T'^ + 7 E Pn^nixV + ^ E u^{m^rr\2^p'y] + . . . | 



(20) 



where the ellipsis represents an infinite string of commutators involving the and M matrix. The commutators 
will always contain at least one power of M making them 0{e) or higher. Each of these matrices can in turn be 
written as sums of their Fourier coefficients multiplied by their appropriate plane wave functions as done above. Since 

plane waves possess the property that u„-^{x)un2{x) ■ ■ .Un,„{x) — w„j+„2H \-n„Xx) one can expand all the sums and 

gather all those terms which are multiplied by u„j+„2H hri,„ (2;) = w„j. (x) for specific ni, 712, . . . and redefine whatever 

fmiction of the pj^'s and (/f^'s that correspond to it, as for rik 7^ 0. 

Given the change of variables in eqn. (fT6| . one must shift the measure of integration from PS to PE and VU^. 
In doing so, one needs to calculate a Jacobian factor that is derived from a Fadeev-Popov procedure that is utilized 
in ref. [5] or equivalcntly from a determinant of the metric of the manifold the fields live on as in ref. 'j^. After 
performing the field redefinition, this factor can be expanded in powers of the non-zero modes of the pion fields: 

[Vi:]=[VU^] Vt (l + 0{M^)y (21) 

For our calculations only the leading order piece is needed. When including nuclcons the ^(x) field is used to describe 
the mesons in addition to E(a;), such that = E. Under SU(2) left and right transformations, E LYiB) =^ ^ 
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LS^V^ or ^ — > VS^K'' while the nucleons transform according to iV ^ VN such that V is an SU(2) matrix that is a 
function of L, R, ^, and x. Under SU(2) vector transformations, V reduces to a spacetime independent element of 
SU(2)y to ensure the correct transformation property of the nucleons. V is defined implicitly through the following 
relations: 



V = VLY^ B (22) 



= \/LE^ (23) 

In building chirally invariant terms for the Lagrangian one must use both transformation rules as stated above. It is 
for this reason that the change of variables needed to separate the zero modes from ^ are given by the following: 

C UiV^ (24) 
C ^ ViU (25) 



and therefore = uWutU and V = VUUJ t/^l^ where L = R'' =U. 

In order to compute the square root of a matrix, assume that there exists a matrix D such that = UTJJ . Since 
the non-zero modes are suppressed by a power of e, the form of D is: 

2 

D = u+jA-^B + ... (26) 

where the factors of e are made explicit, factors of i/ f were put in for later convenience, and A and B can in general 
depend on the pj^ 's and 0° . One can expand UYiU order by order in e and set up a matrix equation to determine A 
and B: 

2 

^ U'^ + '4{UA + ylJ7) - ^ (UB + BU + 2A^) +... (27) 
./ 

= U^ + —UApT''U-—U(ArT''fU + ... (28) 
/ / 

Hence, one is able to construct D to the desired order and V can be readily obtained. Using this to calculate V 
one obtains 

2 

V = 1 + + AU'') - ^ {aI^ - 2AU^M + BU^^ + 0{e') (29) 

In general A and i? are complicated functions of the parameters of U and M . The form of the matrix A is provided 
in appendix D, while the matrix B is not necessary for the calculation at the order considered. 

Under the field redefinition using eqns. (jl6i I24i I25p a few terms in the Lagrangian will be examined to see how 
they are altered. For the pion-nuclcon coupling term given by g"^^ S^A^N, the redefinition causes VAf^V^ 
and so expanding the first couple of terms: 

g'^^^N^Sf^A'^N g^^'^N^S^ (^-^d'^M + ^ |[a^M, M] + AU^d'^M - ^^'MUA^^ + O (e^)^ N (31) 



It is alarming that Af^ now appears to contain terms which are even in the non-zero mode pion fields since this would 
violate parity. Generically denoting a zero mode pion field as tt^, then by comparing the U matrix in the form of: 



(32) 



U = exp 

with that when it is parameterized using hyperspherical coordinates^: 

U = lcos(V') + i(Tisin(?/')sin(0)cos((/)) + i(T2sin(V')sin(0)sin(0) + i(T3sin(?/')cos(6') (33) 
leads to the relations: 

^ = V27Acos(0) 

-y- = -0 sin(6') cxp(— i(/)) 

—J- — 4' sin(6') cxp(z0). 

Substituting eqns. (|33[ I30p in eqn. (j3ip . Taylor expanding the functional dependence in i/), and then using the 
relationships above yield terms with an odd number of non-zero and zero mode pions multiplied by powers of V'^- 
Because ijP' — (|(7rQ)^ + tt^tt^)//^, this term will always be even under parity. Hence, the entire contribution is 
therefore explicitly even under parity. This has been checked to 0(1//"^) and will hold for all orders in 1//. 

The relationship between what are called the non-zero modes under the field redefinition (E UYiU) versus the 
non-zero modes in E must also be carefully explained. For clarity, an expansion of Af_i to three pion fields with exactly 
one nonzero mode pion will be considered. In the expansion of the pion-nucleon Lagrangian term in eqn. (j3ip the 
0(e) term seems to be devoid of zero mode information. Indeed, from the arguments just given, it is reasonable to 
assume that there should exist some function of the zero modes which accompanies it. However, if we again use the 
Baker-Campbell-Hausdorff formula to relate the original parametrization to the one in terms of the field redefinition 
one finds that (where M' encapsulates the non-zero modes before the change of variables and </) encapsulates the zero 
modes): 

1 



6/2 



(34) 

where M is defined by E = exp ( 2iM/ f] . The definition of M will contain an infinite number of commutators, all of 



A precise definition of the coordinates used can be found in appendix C. Normally there is a radial coordinate which we denote by \b\ 
that must exist to parameterize in terms of hyperspherical coordinates. However as shown in the appendix, this factor will be sent to 
one by the presence of a delta function so we explicitly do not write it for this reason. 
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which will be suppressed by a factor of e or more. With this definition one can find 



1 



^d^M' + (fd,Sl' - 2(t)d^M'(^ + 9pM>2) + . . . . (35) 

This result is identical to the expansion of A^_i prior to making the change of variables. In making the change of 
variables all that happened was that some of the zero mode information was absorbed into the definition of M . 
Therefore it is more appropriate to view M not as strictly the non-zero modes of the pions but as an 0{e) suppressed 
quantity with which one can utilize perturbation theory. 

For the mass term of the pions, the quark masses will be taken to be equal, hence the change for the mass term is: 

xtvm,Tr{e + {£}?) -> A^l/m,Tr(;72e2 + (t/t)2(et)2) 

~ X^-^VmqT^{U^ + {U^f)+0{e') . (36) 

The first term on the RHS of eqn. (j36p counts as an order one term and cannot be expanded. It will therefore have 
the interpretation of a probabilistic weight function when considering group integrals over U^. 

Keeping only those terms of the axial-vector current which will contribute to the order at which we work, under 
the redefinition, it goes to: 

= 9a^NS^V {Ut-U^ + C/V"t/) V^N. (37) 
One other equivalent method of separating the zero modes circumvents the use of the V matrix entirely. With this 



method one simply determines the form of using the same method outlined above and uses this directly 

in the Lagrangian. It has not been utilized in this paper because the redefinition in terms of the V^s allows one to 
immediately write down how terms in the Lagrangian are shifted under the change of variables since it mimics a chiral 
transformation. 

In applying the method described above for integrating over the zero mode portion of the pion fields, care must 
be taken in the use of the V and matrices. Specifically, one must be sure that the expansion in powers of M has 
been taken far enough to include all of the appropriate zero modes in the calculation. A good example of this is the 
calculation of the nucleon mass in the e regime performed in ref. 11 1 . To order the graphs that contribute to the 
nucleon mass are given in fig. ([3]). In fig. EJa), [3]Jb), [3Ke), and[3]Jf), the pion couplings are due to one pion derivative 
interactions from the g^'' and (?a7V Lagrangian terms, so to produce these graphs one needs only V ~ 1. 

The contribution from fig. (l3Ud)) arises from the terms in the higher order Lagrangian in eqn. ([5]) proportional to 
Al and Afj.A'^. Under the change of variables these become Afj.A'^ VA^A'^V^ ~ V{d'^nd^n)V\ where again only 
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FIG. 3: Graphs that contribute to the nucleon mass up to order The square vertex is an insertion of the kinetic energy 
operator while the circle is the direct zero mode mass contribution. 



V — 1 is required. Finally at this order in e is fig. (IHc)), which stems from the Lagrangian term 

XCiNTi[mgT,^ + h.c.]N ACi7VTr[m,t/^E^C/^ + h.c.]N 

= ACi7VTr[mg(C/^)2 + /i.c.]iV+ ... 



(38) 



which explicitly carries zero mode information. This term is evaluated in the same manner as the correction to the 
pion mass and results in|ll| 

X'is) 



TUB ^ rriB-^miCi 



2X{s) 



X{s) 



h{2s) 



(39) 



where s = ^rn^f V and the /„(s)'s are modified Bessel functions. The pion mass also has a shifted value due to 



zero-mode contributions such that mz 



mlis) 



11[, where mj(s) = mlX'is)/2X{s). Therefore all TO^'s appearing 



in the computation of the matrix element must be replaced by to^(s). All of these results are identical to those in 



ref. 



Hi 



ll| and gives the same result. However, if one 



To order the method presented above coincides with that of ref. 
wished to carry out the calculation beyond 0{e'^) important differences will surface. An inspection of the original 
Lagrangian will lead to graphs with multiple zero mode loops emanating from the ■nNN vertex and similar vertices. 
Fig. (|4|) displays an example of a graph that will contribute at O(e^) to the nucleon mass (only one zero mode loop 



111 ] some of this 



ll| will be 



is drawn, though an arbitrary number may exist as they will all be of the same order). In ref. 
information does not exist. While correct to the order they worked, the field redefinition used in ref. 
incorrect beyond 0{e^). The inconsistency resides in the field redefinitions and ^ij. The labels of L and R are not 
to be taken literally, but rather are a device to recall the two different ways ^ can transform in order to build chirally 
invariant Lagrangian terms. What is true is that there is only one field, namely ^, such that = E. This must 
hold true both before and after the field redefinitions. Initially, the E field can be thought of as a generic point on 
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the SU(2) manifold that is a finite distance away from the identity. Defining = E imphes that ^ is located on the 
midway point of the geodesic connecting E to the identity. The field redefinitions of U£,l and S^ji U^£,ji lose 

this interpretation. The only way to sensibly interpret the meaning of V UT,U is to appeal to the expansion given in 
eqn. 1^. 
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FIG. 4: One graph that contributes to the nucleon mass at order e^. An arbitrary number of zero mode loops can exist at the 
N{n^)N vertex. 

Turning to the evaluation of the matrix element of the axial- vector current, the part of the action containing the 
sources is expanded keeping only the linear term in each. The matrix element to be evaluated is then of the form: 

< N,ip)\j;,\mip) >=^ f [VmVN] Nlip) N,{p) e-f""^ ^ (40) 



2o. 

where Zq is the partition function free of sources and b and c denote particular choices of nucleons from N"^ — {p, n) 
in the initial and final states. 

N / N A/ N / \ N 
> > ® > — ^ * * 1 t ® > 

(a) (b) 



N / \ N N / \ N 



=®-^— i — > — — ». i Q t y - 

(c) (d) 



N / \ N N , \ N 



(e) (f) 



N \ ) N 
* >^®^> *■ — 



(9) 

FIG. 5: Graphs that contribute at ©(e^) and ©(e'^). 

The graphs in fig. ([5]) represent the perturbative contribution to the matrix element and each of these graphs has 
a function of the zero modes that accompany it which must be integrated over. As an example of how to integrate 
over the zero modes, consider the tree level contribution: 

^Ai(M-\M(^~^ (")r9A7Q /VV «A /[P^'l Ubmmnc cxp (.ReTr([/^)) 

< N^{p)\j^,\mip) >,.,e = 9^ {2NS,N)iT /[pt/^]exp(.ReTr(t/^)) ^^'^ 

= gfQ{s){2NS^N) (42) 
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1 + S^^jiiy J ■ Here, the N's and N are the spinors satisfying ^ (1 + v'^Jfi) N — N (with a similar 



where Q(s) = ^ (l + 2^^^ 

relationship for N) not to be confused with the nuclcon field. In the s oo limit one recovers the expected answer 
at tree- level, g^^\ 

A contribution which enters at C(e^) is the one-loop contribution from the axial- vector current as shown in fig. 
(ISfg)). One can determine the function of s that will modify the graphs from the zero mode integration. The result 
is: 



Tis) 



/o(2g) 



(43) 



The other contributions at this order are all multiplied by Q{s). The functions Q{s) and T(s) are plotted in fig. 




FIG. 6; Zero mode factors Q{s) and T{s) 



RESULTS 



Putting together all of the contributions from the graphs in figure ([5]) and using the results from appendix A the 
result for the e regime to 0{e^) is: 



r,(.) = gfQis) 



Q{s) 



{o)Tis) ciL 
9 a 



100 



(0) 2 \ ( _ caAL^ 



+ 3(3^) 



47r 



-I- 



Q(s) 27r ■ V243^^^ ^ 3^^ ^^^J^l^^^^ 
Co \ 32 (0) 2 ( ciL C2AL^ co 



Co 
2tob 



(44) 



In achieving this result it has been demonstrated in appendix A that the finite time direction corrections in the e 
regime do not enter until 0{e^), thus significantly reducing the computational burden when working in the e regime 
to orders lower than this. 

The result for the e' regime utilizes the graphs from fig. ^ as well as from fig. ([7]). The graphs in fig. (O are 
adapted to the e' regime by letting s 00, which gives Q{s) — + 1 and T(s) 1, and adding back the diagrams with 
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f \ W N / 

A/ V ' A/ » > '%s * *^ 

— » • 

(a) (b) 



N / \ N N / \ 



i — » — — * > jga > ^ 



(c) (d) 

N 4 ^ N 
* >-<S^ *■ 

(e) 

FIG. 7: Additional e' graph contributions. Figure (a) enters at 0{e'). Figure (b) at 0{e''^). Figures (c)-(e) are representative 
of the ways that a zero mode loop can be added to the graphs in fig. [5] in order to create 0{e' ) contributions. 



spatial zero modes: 
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where the subscripts refer to the order in e' considered. In both the e and e' results given above expansions in m.,jL /2t: 



\ An 



(45) 



and AL/27r have been taken as these quantities count as order e and e' respectively. 



VI. DISCUSSION AND CONCLUSIONS 



In this paper we have discussed the various calculational regimes that are relevant to current lattice calculations. 
In fig. ^ the approximate extent of these regimes have been demonstrated for different volumes and pion masses, as 
calculated using the tadpole diagram as a guide. However, it was stressed that these outlines are not to be interpreted 
as strict boundaries, but rather as midpoints in the smooth transition between the different regimes. From the figure 
it is clear that as one moves toward a more hypercubical lattice (with the temporal direction the same length as the 
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spatial extent of the volume) at smaller volumes and more physical pion masses that the e and e' regimes will be the 
most important. 

Ultimately one would like to use the results from the e and e' regime in eqns. (HH) and to fit to lattice data for 
gA and extrapolate to the physical pion mass value. The majority of the data that exists on lattice calculations of gA 



can be found in refs. 
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20|, 



21 



22 



23 



24| . Most of these calculations of the axial current have been performed in the 
e', (5, or p regimes according to the separation given above. While there is a significant amount of data available from 
these calculations, in order to avoid working in the S regime as defined above, only results from volumes significantly 
greater than (2 fm)^ should be used. However, the number of points that strictly satisfy this condition, as well as the 
condition from eqn. (jl2p that m.^L/2'K < 1, is too few for any reasonable fit to be performed. 

In ref. i2^] a fit to the LHPC Collaboration lattice data was performed using a p regime calculation of gA, as 
most of these points fall within that regime. The most effective way to efficiently reduce the statistical error in 
the extrapolation of gA at the physical point (where one wishes to compare to experiment) is to perform lattice 
calculations at lower pion mass (< 350 MeV). Generally speaking, one low mass calculation may easily be equivalent 
to three or four high mass calculations in terms of ability to reduce the statistical error of the extrapolation. While 
performing lower pion mass lattice calculations is often computationally expensive, increases in the speed and number 
of computers available will soon put these low mass calculations within reach. In addition, a high pion mass moves 
one further away from the regime of ChPT, adding to statistical considerations concerns about the applicability of the 
effective field theory, regardless of the regime one is working in. Clearly, the way forward involves lattice calculations 
at lower pion mass, where the e and e' regimes will become increasingly important. Our results allow coherent analysis 
of such data. 

In the e regime the technique of collective variables was used to separate the spacetime independent zero modes 
from the spacetime dependent non-zero modes. Before the split into zero and non-zero modes, one is able to see that 
in graphs where non-zero modes are connected to external states one can have arbitrary numbers of zero mode loops 
present. These contributions are not lost in the field redefinition. To usefully separate out the order one components 
of a given contribution it was demonstrated that one must perform the change of variables that is defined by S = UYJJ 
as in eqn. (|16p . To correctly include baryon interactions after this change of variables the ^ objects must be treated 
carefully with the change of variables given in eqns. and ([^5)1 . We have described a method for doing this that 
differs from previous methods 11 1 . 

One feature of the method proposed for dealing with the 0(1) contributions of the pion fields is that the integration 
over the zero modes must be performed on a diagram by diagram basis. It may seem legitimate to integrate out the 
zero mode contributions and return a Lagrangian where the relevant coupling constants have been replaced by the 
appropriate functions of s, that have been calculated "a posteriori," and with all insertions of the pion fields containing 
purely non-zero mode information. Doing this, however, would be incorrect. Evaluating higher order diagrams will 
lead to different zero mode functions multiplying their respective non-zero mode Feynman diagrams. An example of 
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such a diagram is given in fig. ([8]). These zero mode functions, in general, will not be encoded in a coupling constant 
redefinition; it would be misleading to label gA as gA{s)- As a result, one cannot write down a local Lagrangian after 
the zero modes have been integrated out. To perform the actual integrations on low order diagrams with relatively 
simple zero mode contributions one can cither use hyperspherical coordinates or use the symmetry properties of the 
SU(2) group to obtain an answer. However, with more difficult and complicated higher order diagrams hyperspherical 
coordinates are the most effective tool, as shown in appendix C. 



/ 

N I 



-0= 



FIG. 8: A diagram for wliicli an integration over the zero modes does not correspond to a simple redefinition of coupling 
constants. 

To conclude, the matrix element of the nucleon axial-vector current has been calculated within the e and e' regimes 
to 0{e^) and 0{e'^) respectively. The e' regime result agrees to NLO with previous calculations |l^. A method 
was put forward to calculate the nonperturbative contributions of pion zero momentum modes in the e regime in the 
presence of baryons. This method can be used in any calculation, diagram by diagram, at all orders in the e expansion. 
In calculating gA it was discovered that the temporal direction could be approximated as infinite, with the errors 
introduced entering only at the Oie^) level, as shown in appendix A. This significantly reduces the computational 
burden of calculations in the e regime. 
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APPENDIX 

APPENDIX A: INTEGRAL & SUM MANIPULATIONS 

The integrals used in the calculation of the results in Euclidean space are: 

/rld„ 1 
7-^7^ ^ o (A-1) 

/d'^q \q\'^ 1 

{2n)'^ [tqa - A + q'^/2mBY q'^ + mj^ 
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/cfiq 1 1 |q|^ 

In the e' regime one normally integrates over the go component and forms a sum over the remaining finite volume 
spatial components, while for the e regime a sum over all four components would be formed. However, the e' regime 
can be recovered as a limit of the e finite volume manipulations and so one can focus on those. 

The sums that arise during our calculation are defined within a dimensional regularization framework so that these 
should merge smoothly with the infinite volume pieces as L ~* oo. Toward this end the sums must be defined by the 
relation 

LsA. J (2^)3+^ J (2^)d-i ^^^^ 

where the cutoff dependence of the sum is removed by that of the integral. With this definition, the dimensionally 
regularized integrals will be the infinite volume pieces of the diagram while the remaining sums are the finite volume 
corrections. The DR pieces will all be proportional to or and so will technically enter at order or e'^. 
However, in order to smoothly merge with the infinite volume result these higher order terms must be included along 
with a fi dependent counterterm. These DR pieces are calculated in the normal way 16 1. 



For the temporal sums found in the e regime, the Abel-Plana formula is needed and is given as: 

lv-./27rn\ r dz ^, ^ / f[z) \, / f{z) \, , ^ 

^l^f [ — ) = ^/(^) - '^""^ [ e^l3z _ 1 ) llowcrplane +»RcS [—-pT^ j luppcrplanc (A-5) 

which is valid if f(z) has no poles on the real axis. The first term on the right hand side is the infinite time direction 
piece used in e' calculations. One can thus look at eqn. (|A-5[) as the infinite time e' piece plus finite time e corrections. 
Using eqn. (lA-SP the integral from eqn. (|A-ip becomes: 

d^q 1 



= —y 

fir 3 



(27r)4 g2 +p + m2 

1 



/3i3 ^ (27rno//3)2 + (27rn/L)2 + ml 

(27rno//3)2 + (27rn/L)2 + TO2) ^^"^^ 

Now let P — aL for some value of a where one expects 1 < a < oo. It is now possible to expand the above expression 
in powers of m^L as mT^L ~ with a > 1. Splitting off the zero mode part of this expression and expanding in m^L 
gives 

1 1 



i?o(m.,/3,L) 2 ^^^3 



7.3 ^ 4tt \n\ L3 ^ V 27r Inl e^^"^!"! - 1 

where the infinite time direction term is the first term on the right and the finite time correction is the second term. 
Looking at the sums present in each term, the finite time correction will converge very quickly due to the presence of 
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the exponential. In fact, even for just the first term in the sum the finite time piece will be suppressed by a factor 
of 267 under the most favorable conditions (a = 1). Given that current lattice calculations return e > 1/4, the sum 
for the finite time corrections of the one pion loop is effectively at or higher, and so for the single pion loop, the 
finite time corrections are effectively at and can be safely ignored. Thus: 

i?o(m.,/?,L) = 1-L_ (A-7) 



^ 47r In 

for both the e and e' regimes. 

For eqns. (|A-2p and (|A-3|) the same procedure is used that was used on eqn. (|A-1|) with the only modifications 

being that: (1) the spatial zero modes vanish due to the integrand's proportionality to and (2) the presence of 

the baryon mass in eqn. (|10p must be used in computing the finite volume residues. The pole is in the lower half 
plane and will give a finite time correction of 

- iRes ( -ly^ ^ ) I lower oc „,.y , « (A-9) 



as tob ~ A 



ll| . The nucleon sector will give the same result with A ^ 0. Because of this fact which arose due 
to the finite volume boundary conditions on the relativistic field from eqn. (jlOp . the finite time corrections due to 
the nucleon and decuplet poles do not contribute. From the above results, the finite time corrections from the pion 
propagator poles do not contribute at the order we work. Therefore, to order 0{e^), the finite time contributions to 
the sums are negligible and only the infinite time direction contributions are needed. 
Calculating eqn. (jA-2p explicitly and expanding in m^L and AL: 



J(A,to^,/3,L,tob) = 7771 X! 



(27rn/i)2 1 



/3L3 ^ (i(27rno//3) - A + (27rn/L)2/2mB)2 (27mo//3)2 + (27rn/L)2 + ml 
I ^ ( L I AL2 1 1 \ 4, 



For eqn. (jA-3P it is easier to generalize the integral: 



iV'(A,B,m,,/?,L,mB) - / ^4 , r , ^ ^ //o n l 2//0 V (^-H 

J (27r)* 9^ + 771^ - ^ + q'^/i2mB) iqo - B + q"^ / [2mB) 

N{A,m^,P,L,mB) = iV'(A, 0, m^, /3, L, ms) (A-12) 



Using similar methods to those used to evaluate the function J and taking A A and B gives 

For each of the resulting equations one is left with a spatial sum over various powers of \n\. Several of these sums 
contain divergent pieces which must be removed and placed in appropriate countertcrms. Defining these sums as the 
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limit of an analytically continued sum from non-integer powers of \n\, one can indeed obtain the finite pieces that are 
necessary. This results in 



31 



3: 



C2 



Cl 



Co 



y— 
y- 



-8.913633 



= -2.8372974 



With these coefficients the sums above become 

Ro{m^,l3,L) 
RiiimTr, 13, L) 
J{A,mTr,l3,L,mB) 
N{A,m-^,P,L,mB) 



1 1 



2m^L3 
1 ciL 

I3l7 



C2AL2 




47r2 


2mB 


C2AL2 


h '° 


87r2 


2mB 



(A-14) 

(A-15) 
(A-16) 
(A-17) 

(A-18) 



APPENDIX B: CONTRIBUTING DIAGRAMS 
1. Order e' Graphs 

For the graph in fig. [Zla) the zero next to the pion loop indicates that we are only interested in the zero momentum 
mode of this pion loop. This graph contributes to 0{e'). Its contribution is 



20^°^ 
/ 



(B-1) 



2. Order and e' Graphs 



The graphs in fig. [3a)-(d) yield 



ieldQ: 



2b 



^2c = ^2d = 



100 gAA9ANQis) 

243 p 

32 gfgi^Qis 



J{A,m^,P,L,mB)NS''N 
2 N{A,m^,(3,L,mB)NS''N 



27 p 

The two graphs which contribute to field rcnormalization arc given in figs. [UJe) and (f). 



i 2p. — — - 



2/ 



2 p 

3 P 



J{0,m^,(3,L,mB)NS''N 
JiA,m^,P,L,mB)NS^'N 



(B-2) 
(B-3) 
(B-4) 

(B-5) 
(B-6) 
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Figure [5jg) gives: 



and fig. [7{b) give a contribution of: 



^i^P^i?0(m,,/5,L)iV5'^7V (B-7) 



(0) 

^[R^{m^,l3,L)YNS^'N (B-8) 
where the function Rq){mT^, (3, L) is the non-zero mode portion of the function Rimj^). 

3. Order and e'^ Graphs 

At order e'^, there exist contributions from graphs such as those shown in fig. [Tljc) and (d). To generate these 
graphs one adds a spatial zero mode loop to each vertex in fig. [H^a). As these contain spatial zero mode pion loops 
they contribute only to 0{e'^) and not to 0{e^). Adding a zero mode pion loop to the 0{e'^) graphs at every vertex 
on every graph, including the field renormalization graphs, will generate the majority of the 0{e'^) contributions. 
The effect of adding these pion zero mode loops will be to change the coefficient due to that vertex and multiply the 
previous expressions by the pion loop integral i?o("^7r)- None of the Clebsch-Gordan coefficients or the spin operator 
contractions will change from those used at 0{e'^). To find the change to the vertex coefficients one needs to look at 
each vertex in turn. 

By carefully expanding the vertex terms in the Lagrangian to three pions one can derive the change that occurs to 
the graphs found when a zero mode pion loop is added to figs. (lHa)-(f)). Upon completing the expansion to three 
pions, such a change simply multiplies the graphs from 0{e'^) by a factor of —■^RQ{mTr, (3, L). In addition, one can 
add such a zero mode loop to either end of the graph so that there are two distinct configurations. Hence, adding a 
zero mode pion loop to the Lagrangian vertex will give an 0{e'^) contribution that looks like 



100 2 

^.gATVffAA J (A, TO^, L, mB) 



32 



NS^'N (B-9) 



0(e'^) dia grams are also produced by adding a zero mode pion loop to the vertex due to the current for each of 
the graphs in fig. (O as well as the contribution from fig. El^b). This includes adding a zero mode loop to the current 
vertex that is implicitly part of the field renormalization graphs, creating an 0{e'^) contribution that effectively 
combines the contribution from fig. EJa) with the field renormalizations of figs. [HJe) and (f). Figures [SJa)-(f) are 
multipfied by a factor of —-pRoiniTr, (3, L), while fig. \SIe) is multipfied by —^Ro{mTr,P,L), fig. [T^b) is multiplied 
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by — 5j2-i?o(?Ti7r, /?, i)- Putting all of these current insertion loops together yields a contribution of 



13 , 9a 



f 



1 



100 



32 



(B-10) 



The final 0(e"^) graphs are those as in fig. [Tlje), where the vertex denoted by the dot is from an expansion of the 
pion kinetic energy term. This will generate a four pion vertex, and by choosing the correct contractions of the pion 
fields a zero mode loop can be added as shown to each of the graphs in fig. ([5]). This gives a contribution of 



32 d X6 



H9ff 



10 



J (0, m^, /3, L, tob) + m^^— J J(0, m,r, /3, L, tob) 
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— gAA9AN + -^9a 9aN 



J (A, m^,f3, L, tob) + ml-—r-J{A, m^,(3, L, ms) 

ami 



» (0) 2 
'g5A Sajv 



7V(A, m^, /3, L, tob) + m^— ^iV(A, m^, /?, L, ttib) 

ami 



NS^'N 



(B-11) 



From an operator point of view, both the e and e' regimes have graphs of 0{e^) and ©(e'"^) that look like that in 
fig. ([9]). The square vertex is an insertion of the operator from the 0{l/mB) Lagrangian terms found in cqn. ([5|l. 
However, these graphs have been automatically accounted for by the choice of propagator when the ©(I/tob) terms 
were included. 



N ^ \ N 

* ' m ® ^ — * — 

FIG. 9: Example of operator insertions from the C'(l/m_g) Lagrangian. This graph would contribute at 0{e''^) and 0{e'^) but 
is already counted by the choice of propagator. 



APPENDIX C: HAAR MEASURE 



The integration measure to be used for the zero mode is [DU'^] (parameterized with hypersphcrical coordinates): 



1 



d'^a 5(a^ - 1) 



(C-1) 



such that U"^ = oqI + id ■ a and the integral is normalized to one. For the calculations in this paper U = bgl + ib ■ a 
is parameterized in terms of hypersphcrical coordinates: |5| G [0,oo], tp G [0,7r], 9 £ [0,7r], and (f) € [0, 27r] giving 



U — |6|cos(V')l + i|6|sin(7/')sin(6')cos((/))(Ti + i|6|sin(V')sin(0)sin(0)cr2 + i|6|sin(V')cos(6')cr3 



(C-2) 
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In order to express the integration measure in terms of the radial and angular coordinates describing U we merely 
have to calculate a Jacobian factor. The result obtained is: 

J [VU^] ^ j d% (16fe4cos2(^)) _ 1) 

= ^ j dhb^ (5(fo2 -I) j d^} de d(t> sin2(V')cos2(V')sin(6') (C-3) 

where in the first line the prefactor of 1 /4 is present so the integral is properly normalized and the parenthetical factor 
is the Jacobian. 

The integral in the numerator of eqn. (|4ip can be done by parameterizing the zero mode variable, U, in terms of 
hyperspherical coordinates or in this particular example it can be done directly. To see this, the integrand and measure 
in eqn. (|4ip must be invariant under multiplication on the left and right by arbitrary constant SU(2) matrices. This 
implies that the integral must have the form: 

W^] [/fa,„(C/t)„, exp(sReTr(f/2)) = A{s)SbcS,nn + B{s)Sb,nSnc (C-4) 



Contracting indices and using the relation J[VU^] exp (sReTr([/2)) = X(s)[J,l33|, and 2 det [/ = (Tr U)^ - Tr C/^ 
(valid for SU(2) matrices) one will obtain: 

A{s) = i (2X(s) - X'(.s)) (C-5) 
B{s) = i(X(5) + X'(.)) (C-6) 
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APPENDIX D: MATRIX DEFINITIONS 

The form of the matrix A defined earher in the paper has the following structure (using the hyperspherical 
parametrization fi'om eqn. 



All = -^^ocos('(/') + -^^ocos(0)sin(?/;) + **'n-_sin(6')sin(-0) + **'n-+sin(0)sin(V') 
v2 v2 2 2 

-^e-''^7r_cos(6l)sin(6l)sin(V')tan(V') - ^e"**7r+cos(6l)sin(6l)sin(V')tan(V') 

+ ^7rosin2(0)sin('(/')tan(V') (D-1) 
v2 



A 



12 



A21 = 



1 3 1 

7r+cos('(/)) - -e^^*'^7r_sin(?/')tan(?/') + -7r+sin('(/')tan(V^) + -e"^**7r_cos(20)sin(V')tan(V') 

+ i^-+cos(26')sin(V')tan('(/') ^e"*'^7rosin(26')sin('0)tan(?/') (D-2) 

4 2v 2 

7r_cos(?/') + 7r_cos^(6')sin(?/')tan('(/;) ^e*"^7roCos(0)sin(0)sin(?/;)tan(-0) 

v2 

+ i^_sin2(6l)sin(V')tan(i/') - ie2^'^^+sin2(6')sin(V')tan(V') (D-3) 

— ^7roCos(-0) + ^7roCOs(0)sin(V') + ^e^**7r_sin(0)sin(V') + ^e*"*7r+sin(6')sin(-0) 
v2 v2 2 2 

+ ^e~"^7r_cos(6l)sin(6')sin(?/')tan(V') + ^e*"*7r+cos(6')sin(6')sin('0)tan(V') 

— ^7rosin^(6')sin(V')tan(V') (D-4) 
v2 



[1] E. Jenkins and A. V. Manohar (1991), talk presented at the Workshop on Effective Field Theories of the Standard Model, 

Dobogoko, Hungary, Aug 1991. 
[2] M. Luscher, Commun. Math. Phys. 104, 177 (1986). 
[3] S. R. Beane, Phys. Rev. D70, 034507 (2004), hep-lat/0403015. 
[4] J. Gasser and H. Leutwyler, Phys. Lett. B188, 477 (1987). 
[5] F. C. Hansen, Nucl. Phys. B345, 685 (1990). 
[6] F. C. Hansen and H. Leutwyler, Nucl. Phys. B350, 201 (1991). 
[7] F. C. Hansen (1990), bUTP-90-42-BERN. 

[8] P. Hasenfratz and H. Leutwyler, Nucl. Phys. B343, 241 (1990). 
[9] A. Hasenhatz et al., Z. Phys. C46, 257 (1990). 
[10] H. Leutwyler and A. Smilga, Phys. Rev. D46, 5607 (1992). 

[11] P. F. Bedaque, H. W. Griesshammer, and G. Rupak, Phys. Rev. D71, 054015 (2005), hep-lat/0407009. 
[12] W. Detmold and M. J. Savage, Phys. Lett. B599, 32 (2004), hep-lat/0407008. 
[13] E. Jenkins and A. V. Manohar, Phys. Lett. B259, 353 (1991). 

[14] S.-L. Zhu, S. Pugha, and M. J. Ramsey-Musolf, Phys. Rev. D63, 034002 (2001), hep-ph/0009159. 

[15] T. R. Hemmert, M. Procura, and W. Weise, Phys. Rev. D68, 075009 (2003), hep-lat/0303002. 

[16] S. R. Beane and M. J. Savage, Phys. Rev. D70, 074029 (2004), hep-ph/0404131. 

[17] A. A. Khan et al., Phys. Rev. D74, 094508 (2006), hep-lat/0603028. 

[18] V. Bernard and U.-G. Meissner, Phys. Lett. B639, 278 (2006), hep-lat/0605010. 

[19] M. Procura, B. U. Musch, T. R. Hemmert, and W. Weise, Phys. Rev. D75, 014503 (2007), hep-lat/0610105. 

[20] S. Capitani et al., Nucl. Phys. Proc. Suppl. 79, 548 (1999), hep-ph/9905573. 

[21] D. Dolgov et al. (LHPC), Phys. Rev. D66, 034506 (2002), hep-lat/0201021. 

[22] S. Ohta and K. Orginos (RBCK), Nucl. Phys. Proc. Suppl. 140, 396 (2005), hep-lat/0411008. 

[23] A. A. Khan et al., Nucl. Phys. Proc. Suppl. 140, 408 (2005), hep-lat/0409161. 

[24] R. G. Edwards et al. (LHPC), Phys. Rev. Lett. 96, 052001 (2006), hep-lat/0510062. 

[25] E. Jenkins and A. V. Manohar, Phys. Lett. B255, 558 (1991). 

[26] M. E. Luke and A. V. Manohar, Phys. Lett. B286, 348 (1992), hep-ph/9205228. 

[27] V. Bernard, N. Kaiser, J. Kambor, and U. G. Meissner, Nucl. Phys. B388, 315 (1992). 



24 



[28] A. F. Falk, B. Grinstein, and M. E. Luke, Nucl. Phys. B357. 185 (1991). 

[29] A. Walkcr-Loud, Nucl. Phys. A747, 476 (2005), hcp-lat/0405007. 

[30] H. Lcutwylcr, Phys. Lett. B189, 197 (1987). 

[31] A. Edery, J. Phys. A39, 685 (2006), math-ph/05 10056. 

[32] M. Luscher, Commun. Math. Phys. 105, 153 (1986). 

[33] M. Creutz (1983), Cambridge, Uk: Univ. Pr. (1983) 169 P. (Cambridge Monographs On Mathematical Physics). 



